eeo^'Sovav  ■ 3iij 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
LINCOLN  LABORATORY 


SOME  RESULTS  ON  PHASE-CODED  WAVEFORMS 
WITH  APPLICATION  TO  REDUCTION 
OF  INTERFERENCE  BETWEEN  RADARS 


D.  H.  PRUSLIN 
Group  32 


1 

1 

i 

\ 

i 

I 


TECHNICAL  NOTE  1978-7 

28  FEBRUARY  1978 


Approved  for  public  release;  distribution  unlimited. 


I 


LEXINGTON 


MASSACHUSETTS 


properties  of  binary  phase-coded  radar  waveforms  are  presented 
The  results  are  applied  to  a consideration  of  RFI  reduction 


I Mcqsin^ 
iNm 

' DDC 

U/Wryv^iijiypr. 

^^^STIrlCATlOiV 


CONTENTS 


m 

a 

& 

I 


ABSTRACT 


ill 


i 

1 


I . INTRODUCTION 


II.  SOME  EXAMPLES  EXHIBITING  PROPERTIES 
OF  PHASE  CODES 

A.  "Equivalence"  Between  Cross- 

Correlation  Functions  and  Auto- 
correlation Sidelobes 


B. 


Code  Behavior  vs  Doppler  and 
Code  Length 


C.  Interaction  Between  Phase  Codes 
and  Amplitude  Tapers 


3 

3 


12 

18 


III.  APPLICATION  TO  THE  NOMINAL  SITUATION 
DESCRIBED  IN  SECTION  I 

IV.  FINDING  THE  DESIRED  CODES 

V.  OTHER  APPLICATIONS 

APPENDIX 

ACKNOWLEDGMENT 

REFERENCES 


27 


31 


] 


33 

40 

45 

46 


J 


J 


1 


t 

B 


I. 


INTRODUCTION 


I Recently  there  has  been  renewed  interest  in  the  proper- 

I 

I ties  of  phase-coded  waveforms,  in  general,  and  the  possibility 

that  they  could  be  used  as  a means  for  RFI  reduction,  in  particu- 
lar. This  report  uses  the  latter  as  a focus.  However,  much  of 

f 

! • the  information  developed  is  relevant  to  the  consideration  of 

i 

{ other  applications,  and  some  comments  are  made  in  Sec.  V below. 

\ 

The  nominal  situation  to  keep  in  mind  as  a focus  is  the 
I following.  A set  of  radars  operate  multistatically  within  earshot 

of  each  other  in  the  same  frequency  band.  They  employ  weighted 
burst  waveforms  to  achieve  doppler  resolution  as  well  as  suppres- 
sion against  clutter.  Each  subpulse  of  the  burst  is  a wideband 
waveform  for  reasons  such  as  maximizing  discrimination  poten- 
tial, improving  ECM  immunity,  etc.  Systems  operating  at  L-band 
and  employing  150  MHz  subpulses  can  be  thought  of  as  representa- 
tive. 

In  trying  to  say  something  about  how  phase-coding  could 
be  used  to  ease  the  RFI  problem  implicit  in  the  cdjove  situation, 
questions  come  to  mind  in  the  following  areas: 

1.  Should  the  phase  coding  be  applied  as  a modulation  within 
each  subpulse,  as  a phase  taper  across  the  burst,  or  both? 

In  the  first  instance  the  phase  code  becomes  the  means  of 
pulse  compression  as  well  as  a possible  means  of  RFI  reduc- 
tion, and  thus  autocorrelation  as  well  as  cross-correlation 
considerations  are  important.  In  the  second  instance  the 
phase  taper  will  interact  with  the  aunplitude  taper  being 
used.  If  coding  can  be  applied  at  both  levels,  additional 
RFI  suppression  could  be  possible.  In  all  cases,  the  key 
issue  is  the  waveform  behavior  vs.  doppler  shift. 
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2. 


Two  distinct  approaches  to  achieving  suppression  between 
radars  are  possible  in  the  following  sense.  The  first  is 
for  each  radar  transmitter  to  employ  a distinct  code. 

Codes  are  found  by  some  means  such  that  the  cross- 
correlations satisfy  some  criterion  on  peak  or  average 
level  (average  over  range  extent) . The  second  is  for 
some  of  the  radars  to  employ  a common  code.  The  code 
selected  has  a "deep"  suppression  band  over  some  partial 
segment  of  its  range  extent.  By  scheduling  transmissions 
among  this  set  of  radars,  cross-talk  is  suppressed 
by  the  attenuation  of  the  trough  during  the  listening 
window  of  a given  radar.  The  potential  benefit  is  greater 
suppression.  A disadvantage  is  that  if  the  phase  code 
is  applied  at  the  subpulse  level,  the  window  will  be  very 
short. 

3.  Is  there  knowledge  of  how  to  find  sets  of  codes  satisfy- 
ing given  auto-  and  cross-correlation  constraints  or  how 
to  bound  how  many  exist?  Are  there  design  techniques 
that  yield  codes  having  deep  attenuation  troughs,  speci- 
fied peak  levels,  etc.? 

In  Sec.  II  some  exeunples  are  presented  parametrically 
that  pertain  to  points  1 and  2 above.  In  Sec.  Ill  we  briefly 
apply  these  results  to  the  nominal  example.  Section  IV  contains 
a few  words  about  Number  3 above.  Section  V contains  some  com- 
ments on  other  applications.  Finally  an  appendix  contains  speci- 
fications of  the  codes  used  in  the  examples.  All  codes  used  are 
binary  (O-ir)  phase  codes. 
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II.  SOME  EXAMPLES  EXHIBITING  PROPERTIES  OF  PHASE  CODES 


In  the  examples  to  follow,  range  cuts  are  plotted  vs 
"chips,"  a normalized  time  unit.  If  the  code  is  thought  of  as 
an  intra-subpulse  modulation,  this  time  unit  is  the  compressed 
pulsewidth.  If  the  code  is  thought  of  as  a phase  taper  across 
the  burst,  this  time  unit  is  the  subpulse  spacing.  When  quan- 
tities are  plotted  vs  doppler  frequency  a normalized  abscissa 
is  used.  Unity  on  this  axis  corresponds  to  the  reciprocal  of 
the  compressed  pulsewidth  or  of  the  subpulse  spacing,  respec- 
tively, in  the  two  cases.  The  difference  will  be  important  in 
Sec.  III. 

A.  "Equivalence"  Between  Cross-Correlation  Functions 
and  Auto-Correlation  Sidelobes 

Since  less  is  known  about  cross-correlation  properties 
of  phase  codes  than  auto-correlation  properties , most  codes  dis- 
cussed or  tcUbulated  in  the  literature  have  some  specified  auto- 
correlation characteristics . Since  we  are  interested  in  both 
types  of  functions , the  purpose  of  this  section  is  to  examine 
to  what  extent  cross-correlation  functions  behave  more  or  less 
like  sidelobe  patterns  of  auto-correlation  functions.  To  this 
end,  we  will  consider  two  separate  pairs  of  phase  codes.  The 
first  pair  can  be  characterized  as  having  "good"  auto- correlation 
properties , with  acceptance  of  whatever  cross-correlation  prop- 
erties go  along  with  this.  The  second  pair  can  be  characterized 
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as  having  been  in  some  sense  "designed"  to  have  reasonable 
cross-correlation  as  well  as  auto-correlation  properties. 

The  first  pair  is  a pair  of  known 31  chip  codes. 

One  has  the  minimum  peak  auto  sidelobe  of  any  binary  code  of 
this  length.  The  other  has  the  minimum  peak  auto  sidelobe  of 
any  binary  maximal  length  shift  register  code  (see  appendix)  of 
this  length.  Figure  1-a  shows  the  auto  and  cross-correlations 
(zero  doppler  range  cuts).  Figure  1-b  shows  the  peak  and  aver- 
age values  on  a range  cut  of  each  function  vs  the  normalized 
doppler  frequency  of  the  range  cut.  On  the  auto-correlation 
curves  the  mainlobe  or  peak  occurring  at  zero  shift  is  excluded. 
None  of  the  tiiree  functions  happen  to  have  an  extended  trough 
of  high  attenuation.  Hence,  Fig.  1-c  shows  for  each  of  the  three 
functions  the  behavior  vs  doppler  of  a selected  sidelobe  which 
is  low  at  zero  doppler. 

In  Fig.  1-b  the  noticeable  difference  between  the  auto 
and  cross  functions  is  the  degradation  in  peak  and  average  side- 
lobe level  at  the  greater  doppler  offsets  on  the  auto  function, 
and  the  lack  of  degradation  on  the  cross  function.  The  effect 
is  most  pronounced  for  peak  level.  A possible  rationalization 
for  this  difference  is  that  these  two  codes  are  somewhat  atypical 
with  respect  to  auto  sidelobe  level,  having  been  optimized  in 
that  regard  at  zero  doppler.  Thus,  they  pop  up  more  strongly 
off  doppler  while  the  cross-correlation  level,  not  being  spe- 
cially chosen  at  zero  doppler,  is  less  sensitive  to  doppler  offset. 
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This  fragile  explanation  can  be  jarred  by  considering 

roi 

the  second  pair  of  codes.  These  codes ‘ are  101  chips  in  length, 
and  unlike  the  previous  pair,  were  determined  by  a method  which 
atten^ted  to  distribute  energy  evenly  between  the  cross  correla- 
tion level  and  the  auto-correlation  sidelobes.  Figure  2-a  shows 
the  zero  doppler  range  cuts.  Fig.  2-b  shows  the  behavior  of  the 
peak  and  average  sidelobe  levels  vs  doppler  offset,  and  Fig.  2-c 
shows  the  behavior  vs  doppler  of  selected  sidelobes  which  have 
high  attenuation  at  zero  doppler.  The  similarity  between  Figs. 

1- b  and  2-b  is  striking.  The  average  cross-correlation  level  is 
virtually  independent  of  doppler  shift,  while  the  average  auto 
sidelobe  levels  exhibit  some  degradation.  The  peak  cross  corre- 
lation level  degrades  slightly,  while  the  peak  auto  sidelobe 
levels  degrade  more  strongly,  as  well  as  more  strongly  than  their 
average  levels . 

A more  pronounced  effect  can  be  seen  in  Figs  1-c  and 

2- c.  In  each  case  it  is  seen  that  a "low"  isolated  sidelobe 
level  on  the  zero  doppler  range  cut  degrades  strongly  vs  doppler 
con^Jared  with  the  peak  or  average  sidelobe  levels  discussed  above. 
This  behavior  should  be  taken  as  suggestive  of  how  low  level 
troughs  would  behave  (see  next  section) . The  cross-correlation 
behavior  is  very  similar  to  the  auto  behavior. 

The  suggested  conclusions  to  be  drawn  from  the  examples 
of  this  section  (and  others  not  presented)  are; 
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Fig,  2-b,  Peak  and  average  sidelobe  levels  vs  Doppler 
shift  - 101  chip  codes. 


1.  Cross-correlation  functions  experience  less  degradation  vs 
doppler  than  the  modest  amount  experienced  by  auto  correla- 
tion sidelobes,  as  far  as  peak  and  average  levels  are  con- 
cerned. Peak  level  is  the  more  sensitive  of  the  two  quanti- 
ties. It  can  be  speculated  that  this  behavior,  which  seems 
to  apply  when  the  zero  doppler  cross-correlation  level  is 
greater  than  the  auto  sidelobe  level , might  be  reversed  for 
a pair  of  sequences  in  which  the  cross-correlation  level 
were  the  smaller  quantity. 

2.  Greater  degradation  vs  doppler  is  experienced  by  low-level 
sidelobes  (those  corresponding  to  one  or  a few  chips  resi- 
due) on  the  zero  doppler  cuts  of  either  type  of  function. 

The  "typical"  behavior  is  similar  for  either  type,  auto  or 
cross  correlation. 


B.  Code  Behavior  vs  Doppler  and  Code  Length 
With  Sec.  A as  background,  this  section  employs  exeun- 
ples  involving  auto-correlation  functions  to  examine  how  code 
behavior  depends  not  only  on  normalized  doppler  frequency,  but 
code  length.  Just  as  there  is  a big  difference  between  how  real 
doppler  shift  relates  to  normalized  doppler  when  the  phase  code 
is  applied  as  a subpulse  modulation  and  when  it  is  applied  as  a 
burst  weighting,  so  too  there  will  be  differences  in  how  code 
response  varies  vs  normalized  doppler  for  different  code  lengths. 
This  is  so  because  for  a given  chip  width,  a longer  code  provides 
more  time  for  doppler  induced  phase  shifts  to  develop  across  the 
waveform.  The  expectation  is  that  longer  codes  will  degrade  more 
quickly,  and  the  purpose  of  the  excunples  to  follow  is  to  quantify 
this  to  some  degree,  as  well  as  to  demonstrate  some  behavior 
hinted  at  in  Sec.  A. 
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The  three  codes  (see  Appendix)  used  as  exan^les  have  in 
common  the  feature  that  their  zero  doppler  (auto)  range  cuts  have 
an  "extended"  region  adjacent  to  the  auto  correlation  spilce  at 
the  origin,  in  which  the  sidelobe  level  is  "very  low."  Very  low 
means  that  the  voltage  sidelobe  levels  in  the  region  are  small 
multiples  of  1/N,  where  N is  the  number  of  chips  in  the  code.  In 
addition,  the  codes  all  have  "reasonably  low"  overall  sidelobe 
levels,  i.e.,  the  range  cuts  looJc  like  a spike  at  the  origin 
sticking  out  well  above  the  proverbial  "rumpled  rug"  elsewhere. 
The  code  lengths  are  32,  255,  and  1023.  These  values  pretty  much 
span  from  a length  that  would  be  used  across  a burst  (several 
tens)  to  a length  applicable  to  a sub-pulse  modulation  (in  the 
hundreds) . The  32  and  1023  chip  codes  have  not  received  specific 
attention  in  the  open  literature.  The  255  chip  code  is  mentioned 
in  Nathanson's  text  book. 

The  range  cuts  of  these  codes  are  shown  in  Fig.  3-a. 
Figure  3-b  shows  the  overall  sidelobe  behavior  of  the  codes  vs 
doppler  shift,  while  Fig.  3-c  shows  the  behavior  of  the  low 
amplitude  trench  off  doppler.  The  following  observations  can  be 
made  regarding  Figs.  3-b  and  3-c. 

1.  The  suppression  level  achievable  increases  with  code  length, 
both  in  terms  of  overall  peak  or  average  level  and  of  trench 
depth. 

2.  At  zero  doppler,  the  overall  peak  levels  are  near  1/N  in 
power  units  with  average  levels  somewhat  lower.  Far  enough 
off  doppler  the  suppression  degrades  so  that  the  peak  and 
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Autocorrelation  trench  levels,  peak  and  average 


average  levels  bracket  1/N.  Pairs  of  codes  exist  for  which 
both  auto-correlations  and  the  cross-correlation  approximate 
these  levels.  If  the  overall  auto  levels  are  pushed  lower, 
the  cross  level  will  be  higher  and  vice  versa,  l2] 

2 

3.  The  minimum  achievable  non- zero  trench  level  is  1/N  (power 
units)  corresponding  to  a 1 chip  "residue."  The  definition 
of  the  extent  of  a trench  is  arbitrary  , but  typically  residue 
levels  of  several  chips  will  be  included,  which  iSjWhy  some 
of  the  peak  levels  in  Fig.  3-c  are  higher  than  1/N^  (a  4 chip 
residue  appears  12  dB  higher  on  the  plot) . Off  doppler,  the 
degradation  is  of  much  greater  magnitude  and  occurs  much 
sooner  than  the  degradation  in  overall  level  for  a given  code 
length.  The  trench  levels  rise  to  approximately  the  same 
values  as  the  overall  levels. 

4.  The  longer  the  code,  the  sooner  it  degrades  with  doppler. 
Breakpoint  frequencies  scale  about  inversely  with  code 
length. 

In  addition  to  the  sidelobe  properties  discussed  above, 
there  is  another  characteristic  which  must  be  considered  when  the 
autocorrelation  function  is  important,  i.e.  , when  the  code  is 
used  as  a subpulse  modulation  for  pulse  compression.  That  char- 
acteristic is  the  behavior  of  the  auto-correlation  spike  at  zero 
range.  The  zero  doppler  auto  range  cut  looks  like  a spike  riding 
on  a rumpled  rug  in  this  case.  We  have  been  examining  how  the 
rug  behaves  vs  doppler  shift.  The  other  degradation  that  must  be 
checked  for  is  the  disappearance  of  the  zero  range  spike  with 
doppler,  leaving  only  the  rumpled  rug.  This  effect  corresponds 
to  nothing  more  than  the  "basic  doppler  resolution"  of  a rectan- 
gular pulse  N chips  long,  since  the  phase  code  has  no  effect  on 
the  doppler  response  at  zero  range.  This  is  easily  seen  from 
the  definition  of  the  ambiguity  function: 
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X(T,f) 


00 


(1) 


= J U(t)u*(t+T)e"^^^^^dt 


where  u(t)  is  the  complex  signal,  phase  coded  in  our  case.  The 
dependence  of  the  spike  amplitude  on  doppler  is  then 


00 

x(0,f)  = J |u(t)  I V^^'^^^dt 
^•00 


(2) 


i.e.,  the  Fourier  transform  of  the  magnitude  squared  of  the  complex 
signal.  The  phase  code  does  not  affect  this,  and  the  spike  ampli- 
tude varies  as 


lx(o,f) I 


Sin  Nirf 
Nirf 


(3) 


along  the  normalized  doppler  axis.  Figure  3-d  shows  this  depen- 
dence . 

C.  Interaction  Between  Phase  Codes  and  Amplitude 
Tapers 

This  section  contains  examples  pertaining  to  the  com- 
bined use  of  a phase  code  and  an  amplitude  taper.  We  now  have  in 
mind  working  at  the  burst  level.  The  "original"  amplitude  taper 
is  employed  to  obtain  the  desired  burst  doppler  response,  Q- 
function,  etc.,  i.e.,  to  achieve  suppression  against  "slow"  clut- 
ter. The  question  is,  what  will  happen  to  this  doppler  response 
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if  a phase  code  or  taper  is  applied  across  the  burst  to  achieve 
some  suppression  against  signals  received  from  other  radars  trans- 
mitting nearby? 

For  an  example  we  take  a 31  pulse  burst.  Figure  4-a 


shows  the  zero  doppler  range  cuts  for  three  types  of  (conplex) 
weighting.  One  is  the  use  of  Hamming  amplitude  weighting  alone, 
which  results  in  a 30  dB  Q-function  and  better  than  60  dB  suppres- 
sion against  slow  clutter  on  the  first  several  range  ambiguity 
splines.  The  second  is  the  use  of  a phase  code  alone,  in  parti- 
cular the  31  chip  Delong  phase  code  having  minimum  peak  sidelobe, 
previously  referred  to  in  Sec.  II-B.  The  third  is  the  combined 
use  of  these  anplitude  and  phase  tapers. 

The  range  cut  for  the  amplitude  weighted  burst  exhibits 
the  behavior  typical  of  weightings  that  yield  good  suppression 
against  extended  slow  clutter;  a "slow"  falloff  of  the  near-in 
ambiguity  amplitudes.  The  range  response  of  the  phase  coded  wave- 
form drops  abruptly  off  the  origin,  a behavior  associated  with 
poor  suppression  against  slow  clutter.  The  underlying  principle 
here  is  the  relationship  which  holds  between  the  Q-function, 
which  measures  doppler  suppression  for  fully  range-extended  clut- 
ter, and  the  zero-doppler  range  cut;  they  are  Fourier  transforms. 
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Fig.  4-a.  Zero-Doppler  Range  Cuts,  31  pulse  Burst. 
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To  get  good  suppression  at  dopplers  far  from  that  of  the  target, 
i.e.,  a low  Q-function  floor  or  "high  frequency"  content  in  the 
Fourier  transform,  requires  a "gentle"  rate  of  decrease  in  the 
high  energy  portion  of  the  range  cut. 

The  range  cut  for  the  combined  an^litude  and  phase 
tapered  burst  is  dominated  by  the  phase  code , since  the  same  se- 
quence of  overlaps  occurs  between  0-chips  and  ir -chips , the  ampli- 
tude weighting  acting  only  to  reduce  the  cancellation  level  (i.e., 
increase  the  ambiguity  level)  somewhat  near  the  origin.  Thus  we 
expect  to  find  degraded  clutter  suppression  relative  to  that  of 
the  original  amplitude  weighted  burst. 

The  degree  of  degradation  depends  on  the  actual  range 
extents  encountered  which  in  turn  depend  on  system  and  threat 
geometry.  If  fully  extended  clutter  must  be  dealt  with.  Fig.  4-b 
which  shows  the  Q-functions  for  the  amplitude  weighted  and  anpli- 
tude/phase  weighted  bursts  , indicates  that  the  burst  clutter 
suppression  is  completely  lost.  If  a bistatic  geometjry  demands 
clutter  suppression  on  only  the  near-in  ambiguity  splines,  the 
picture  is  not  as  bleak.  Figure  4-c  shows  the  doppler  cuts  for 
the  amplitude  weighted  burst  at  zero  range  and  on  the  first  two 
ambiguities.  More  than  60  dB  of  doppler  suppression  is  achieved 
as  a response  floor.  Figure  4-d  shows  the  corresponding  responses 
for  the  amplitude/phase  tapered  burst.  The  doppler  response  at 
zero  range  is  undegraded.  Equation  2 is  the  reason.  The  responses 
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Fig.  4-b.  Q-Functions  vs  Doppler  shift,  31  pulse  burst 


Fig.  4- 
splines 


c.  Doppler  response  on  first  three  ambiguity 
- Hamming  eunplitude  weighting. 
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Pig.  4-d.  Doppler  response  on  first  three  ambiguity 
splines  - combined  amplitude  weighting  and  phase  code. 
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on  the  next  two  ambiguities  are  greatly  degraded  but  still  pro- 
vide 10  to  15  dB  of  clutter  suppression  on  the  average.  It  be- 
comes a question  of  how  much  clutter  suppression  is  needed. 

The  phase  code  used  to  make  the  points  above  was  chosen 
for  no  better  reason  than  that  we  had  it  laying  around,  really. 
While  the  results  obtained  with  it  above  are  typical  of  what 
would  happen  with  most  known  phase  codes,  the  fact  is  that  most 
of  these  codes  are  wrong  for  the  application  of  this  section. 

Most  of  these  codes  have  "low"  autocorrelation  sidelobes.  This 
will  lead  to  a high  degradation  of  the  burst  doppler  response. 

Low  range  sidelobes,  ambiguities  in  this  approach,  are  not  a re- 
quirement. The  original  amplitude  weighted  burst  has  high  range 
sidelobes  anyway.  What  is  needed  is  a set  of  codes  with  high 
auto  sidelobes  and  low  cross-correlation.  This  combination 
will  maximize  RFI  suppression  and  minimize  degradation  of  clutter 

suppression.  Such  codes  are  the  antithesis  of  the  type  generally 

[21 

searched  for  and  tabulated.  There  is  evidence'  that  looking 
for  a pair  of  codes  with  high  auto  sidelobes  and  low  cross - 
correlation  is  not  fighting  nature.  If  such  a combination  were 
found  the  limiting  factor  on  applying  the  phase  code  at  the 
burst  level  could  become  the  relatively  faster  degradation  of 
cross-correlation  level  with  doppler  shift  because  of  the  long 
chip  size  (equal  to  the  burst  subpulse  spacing)  of  the  code. 
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III • APPLICATION  TO  THE  NOMINAL  SITUATION  DESCRIBED  IN  SECTION  I 

As  a nominal  example,  consider  L-band  radars  at  1000  MHz, 
a 50  us  burst  length,  and  150  MHz  bandwidth.  Let  the  ambiguous 
velocity  be  10  knv/sec  to  keep  the  low  velocity  clutter  out  of  the 
mainlobe  of  the  first  ambiguous  doppler  response.  The  ambiguous 
doppler  frequency  is ; 

2v 

^AMB  ^ ^ =66.7  KHz  (5) 

X 

The  corresponding  subpulse  spacing  is 

A = = 15  us  (6) 

AMB 

This  gives  us  3 pulses  in  50  us  which  isn't  much  of  a burst. 
Therefore,  let's  shrink  the  subpulse  spacing  to  5 us,  correspond- 
ing to  an  U-pulse  burst,  and  moving  the  doppler  ambiguity  out 
to  30  km/sec  from  the  primary  doppler  channel.  The  nominal 
doppler  resolution  of  the  burst  is  1/50  usee  = 20  KHz,  or  3 kitv'sec. 

As  for  the  150  MHz  bandwidth,  this  corresponds  to  com- 
pressing the  5 us  subpulse  spacing  to  6.7  ns,  a compression  ratio 
of  750.  The  doppler  ambiguity  or  normalized  doppler  of  unity 
point  corresponding  to  6.7  ns  is  22,500  km/sec,  or  750  times  the 
30  km/sec  value  corresponding  to  the  5 usee  subpulse  spacing  of 
the  burst. 

Now  consider  using  phase  coding  within  the  subpulses 
for  the  joint  purposes  of  pulse  conpression  and  isolation  from 
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undesired  returns  due  to  other  transmitters.  The  phase  codes  are 
then  750  chips  long.  A giant  inferential  leap  from  Figs.  2-b 
and  3-b  indicates  that  a set  of  codes  "can  be  found"  such  that 
with  no  doppler  offsets  involved,  the  peak  auto-correlation 
sidelobe  levels  (pulse  compression  range  sidelobes)  and  peak 
cross-correlation  levels  (RFI  suppression  level)  would  be  1/N 
in  power  units  or  -28.8  dB,  with  average  levels  5 to  10  dB  lower. 
If  the  range  sidelobe  levels  are  too  high,  spectral  weighting 
could  be  applied  in  some  form.  We  have  not  examined  the  effect 
of  this  on  the  cross-correlation  levels. 

As  far  as  doppler  offsets  are  concerned,  one  wants  to 
be  sure  that  the  auto-correlation  doesn't  degrade  over  the  range 
of  velocities  expected  of  legitimate  targets  , and  that  the  cross- 
correlation level  is  maintained  for  RFI  returns  reflected  from 
moving  targets.  7 km/sec  will  bound  the  range  of  velocities 
we  are  dealing  with,  a normalized  doppler  extent  of  .0003.  Again, 
taking  Fig.  3-b  as  representing  what  will  happen,  this  degree  of 
doppler  offset  will  result  in  negligible  degradation  of  peak  and 
average  levels,  although  the  breakpoint  is  being  approached.  As 
for  maintenance  of  the  central  auto-correlation  spike.  Fig.  3-d 
indicates  that  this  is  not  a problem. 

Regarding  the  approach  of  using  a common  phase  code  with 

a deep  trench  and  scheduling  of  transmissions.  Fig.  3-c  indicates 

2 

that  a peak  trench  level  of  1/N  can  be  approached  as  a limit. 
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-57.5  dB  in  this  case.  Allowing  a peak  residue  of  4 chips  in  the 
trench  will  correspond  to  a peak  level  12  dB  above  this  value, 
or  about  -45  dB  with  an  average  level  of  -50  to  -55  dB.  However, 
these  values  will  degrade  more  quickly  with  doppler  than  the 
28.8  dB  value  above.  Using  Fig.  3-c  as  a guide,  we  can  guessti- 
mate a 10-15  dB  degradation,  bringing  the  peak  suppression  corres- 
ponding to  undesired  transmissions  reflected  off  slow  clutter  up 
to  near  -30  dB. 

The  other  problem  is  the  range  extent  of  the  trench. 
Guessing  at  10  percent  of  the  code  length  as  what  could  somehow 
be  achieved,  we  get  75  chips  which  corresponds  to  10  meterslll 
Thus,  in  this  way  of  doing  things  we  have  the  absurd  result  that, 
assuming  all  the  system  timing  could  be  worked  out,  a listening 
window  10  meters  long  would  be  cleared  of  RFI.  This  is  not  very 
useful  for  surveillance. 

What  about  the  second  basic  alternative,  that  of  apply- 
ing the  phase  code  across  the  burst?  The  peak  cross -correlation 
level  would  now  be  only  1/11  = -10.4  dB.  This  could  be  thought  of 
as  an  additional  isolation  achievable  by  applying  phase  codes  at 
both  subpulse  and  burst  levels.  Since  7 km/sec  corresponds  to 
a normalized  doppler  frequency  of  .23  at  5 ps  subpulse  spacing, 
this  value  will  experience  whatever  degradation  occurs  vs  doppler 
for  the  particular  code  used,  if  RFI  is  received  off  of  slow  clut- 
ter. As  we  have  seen,  the  big  problem  is  the  severe  degradation 
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of  burst  doppler  response  on  all  but  the  central  range  ambiguity 

The  trench  approach  would  yield  at  best  a peak  suppression  level 
2 

of  (1/11)  = -20.8  dB,  and  the  extent  of  the  trench  would  now  be 

measured  in  kilometers  rather  than  meters.  An  auto-correlation 
function  of  this  type  would  probably  lead  to  nearly  maximum  deg- 
radation of  the  burst  doppler  response,  and  the  trench  level 
would  degrade  significantly  over  the  range  of  possible  doppler 
shifts  for  RFI  returns. 
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IV.  FINDING  THE  DESIRED  CODES 


As  can  be  noticed  from  Sec.  II  and  Sec.  Ill,  we  have 
attenpted  to  make  our  points  by  referring  to  examples  that  were 
supposed  to  typify  some  kind  of  behavior.  The  reason  for  this 
is  that  there  are  essentially  no  equations  that  relate  an  output 
waveform  characteristic  such  as  a peak  sidelobe  level  to  any 
specification  of  the  phase  code  itself.  Stated  another  way, 
there  are  really  no  design  techniques  that  produce  phase  codes 
with  specified  properties  such  as  auto- correlation  cuts  with 
trenches  of  given  depth,  specified  overall  sidelobe  levels,  etc. 
What  one  finds  basically  are  tabulations  of  codes  (or  genera- 
tional information)  which: 

1.  Were  found  by  enumeration, 

2.  Satisfy  some  condition  on  zero  doppler  auto-correlation 
cut  only,  e.g.,  lowest  peak  sidelobe  for  code  of  given 
length. 

3.  May  fall  into  some  special  subclass,  e.g.  , maximal  length 
shift  register  codes. 

In  terms  of  pairwise  cross-correlation  properties,  not  to  mention 

mutual  cross  correlation  properties  among  codes  in  a larger  set, 

(21 

one  doesn't  even  find  this.  A recent  paper  ‘ ^ describes  an  algo- 
rithm for  starting  with  a pair  of  codes  and  altering  the  relation 
between  the  auto  correlation  sidelobes  and  cross-correlation  levels. 
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The  method  has  some  promise  and  is  an  improvement  over  blind 
search,  but  it  is  still  a long  way  from  a multi-sequence  design 
technique.  Thus  finding  sets  of  codes  with  the  desired  proper- 
ties will  be  a hunt  and  peck  procedure,  and  there  are  an  awful 
lot  of  binary  codes  of  length  750. 

Short  of  a method  for  finding  codes  , it  would  be  use- 
ful to  know  how  many  codes  of  a given  length  and  satisfying 
certain  requirements  on  their  auto  and  cross-correlations  exist, 

or  at  least  bounds  on  this  number.  A lower  bound  has  been  ob- 
f4l 

tamed  on  the  number  of  binary  codes  having  maximum  auto- 
correlation sidelobe  no  greater  than  a specified  value  and 
maximum  cross-correlation  level  between  any  pair  no  greater  than 
another  specified  value.  Unfortunately,  for  values  of  the  para- 
meters of  interest  here,  this  lower  bound  is  zero,  so  no  informa- 
tion on  the  size  of  such  a code  set  is  gained  (the  two  cases  of 
interest  are:  low  auto  sidelobes , say  -30  dB  or  less,  100-1000 
chips,  any  value  of  cross-correlation;  i.e.  , the  use  of  coding 
within  a subpulse; or:  high  auto  sidelobes,  a few  tens  of  pulses, 
any  value  of  cross-correlation;  i.e.,  the  case  of  coding  across 
the  burst) . 
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V.  OTHER  APPLICATIONS 


As  we  have  seen,  phase  coded  waveforms  exist  which 
exhibit  "extra-deep"  suppression  regions  adjacent  to  the  central 
response  on  the  range  axis.  This  observation  has  led  to  some 
conjecture  that  if  the  low  response  level  extended  out  from  the 
range  axis  then  it  might  be  possible  to  encompass  most  of  a clut- 
ter distribution  within  the  suppression  region  by  adjusting  the 
waveform  parameters  and  possibly  by  applying  some  transformation 
to  the  waveform.  One  scenario  of  interest  would  be  designation 
in  a tank  break-up  (TBU)  environment,  in  which  the  clutter  has 
significant  range  and  doppler  extent.  Figure  5-a  shows  qualita- 
tively what  a TBU  clutter  distribution  might  look  like  at  a time 
when  a surveillance  or  bulk  filtering  waveform  would  be  employed. 
Figure  5-b  represents  the  ambiguity  function  of  a hypothetical 
signal  in  which  the  sort  of  response  shown  in  Fig.  3-a  extends 
indefinitely  in  the  doppler  direction,  while  Fig.  5-c  represents 
the  rotation  of  this  response  function  in  the  range-doppler  plane 
as  a result  of  some  transformation  applied  to  the  signal. 

Let  us  deal  first  with  the  rotation  of  the  response 
function.  It  is  conceivable  that  in  some  situations  the  rotated 
form  of  a response  function  could  provide  a better  fit  to  a given 
clutter  distribution.  The  question  here  is  how  must  a phase-coded 
signal  be  modified  to  produce  such  a rotation.  The  answer  comes 
from  the  following  general  relationships. 
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Consider  a complex  radar  signal  x(t),  having  Fourier 


transform 


= J x(t)e“^^’^^^dt 


and  a receiver  matched  to  a signal  w{t)  having  Fourier  transform 
W(f).  The  cross-ambiguity  function  of  this  (generally)  mismatched 
system  is 


Xxw(t,v)  = J*  x(t)w*  (t+T)e~^^’^^*^dt 


Next  consider  a new  radar  signal  y(t)  which  is  the  Fourier 
transform  of  the  original  signal  x(t) ; 


y(t)  = X(t) 


and  a receiver  matched  to  a signal  v(t)  which  is  the  Fourier 
transform  of  the  original  weighting  function  w(t) : 


v(t)  = W(t) 


One  can  then  show  that  the  cross-ambiguity  function  for  the  new 
signal  pair  satisfies 

|Xy^(T,v)  I = IXxv,(v,-T)|  (11) 

• i.e.,  that  the  response  functions  are  related  by  a 90“  rotation 

] 

i in  the  range-doppler  plane. 

To  rotate  the  response  function  of  a phase-coded  wave- 
form and  its  matched  receiver  by  90“  we  must  then  transmit  the 
Fourier  transform  of  this  signal  as  a time  waveform,  and  match 
the  receiver  accordingly.  What  does  the  Fourier  transform  of 
a "typical"  phase-coded  waveform  look  like?  The  answer  in  general 
is  that  it  is  characterized  by  both  a continuously  eimplitude 
modulated  magnitude  function  and  a continuously  varying  phase 
function.  Remember  that  this  magnitude  and  phase  function  would 
have  to  be  transmitted  ^ time  in  order  to  rotate  the  response 
function  of  the  related  phase-coded  signal. 

We  see  that  to  achieve  any  benefit  accruing  to  rotation 
in  a particular  situation  we  must  pay  the  price  of  trading  in  a 
phase-coded  signal  with  its  desirable  properties,  viz,  lack  of 
amplitude  modulation,  ease  of  generation,  simplicity  of  phase 
modulation,  straightforwardness  of  receiver,  on  a new  model  which 
has  the  converse  properties  and  makes  a very  unattractive  candi- 
date for  implementation. 


Let  us  now  return  to  the  question  of  fitting  the  phase- 
coded  response  function  to  the  clutter  distribution.  Assume  that 


we  are  considering  the  same  32  chip  signal,  for  which  the  range  ex- 
tent of  the  trench  is  +6  chips.  To  roughly  encompass  the  range 
extent  of  the  clutter  let  6A  = 6 km,  or  A = 1 km  -►  6.7  ys,  where 
A is  the  chip  width  or  spacing.  Referring  to  Fig.  5-b  which  illus- 
trates the  hypothetical  response  function,  we  observe  that  it  dif- 
fers from  the  actual  response  function  of  this  signal  in  at  least 
two  respects.  Firstly,  the  central  response  peak  at  the  origin 
will  be  replicated  at  zero  relative  range  at  the  Doppler  eunbiguity 
corresponding  to  A,  namely  X/2A,  where  X is  the  radar  wavelength. 

If  this  ambiguity  falls  within  the  clutter  the  desired  target 
response  will  be  swamped,  and  range  cells  to  either  side  of  the 
target  will  be  filled  with  clutter  returns.  Suppose  we  select  a 

j 

radar  frequency  of  2 GHz,  placing  operation  between  L-  and  S-band. 

Then  X = . 15  meter  and  = 11.2  km/sec,  which  is  safely  out  of 

the  clutter.  j 

The  second  discrepancy  between  the  hypothetical  and 
actual  response  functions  is  the  one  which  fundamentally  undercuts  ' | 

any  attempt  at  suppressing  clutter  with  significant  range  and 
doppler  extent  using  the  suppression  region  of  a phase-coded  res- 
ponse function.  This  discrepancy  corresponds  to  the  fact  that, 
far  from  extending  indefinitely  in  the  doppler  dimension  as  in 
Fig.  5-b,  the  suppression  region  has  a doppler  extent  less  than 
1/NA,  the  doppler  resolution  of  the  signal.  Figure  3-c  illustrates 

1 
I 
i 

I : 

i- 
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this.  Thus  rather  than  appearing  as  in  Fig.  5— b,  the  suppression 
region  looks  more  like  that  shown  in  Fig.  5-d.  In  the  present 
example  the  doppler  resolution  is  11.2/32  = .3  km/sec,  so  the 
suppression  region  doesn't  even  reach  the  clutter  in  doppler,  no 
less  encompass  it. 

In  summary,  "special"  properties  of  phase-coded  signals 
can  be  taken  advantage  of  for  clutter  suppression  only  when  dop- 
pler differentials  or  extents  are  qualitatively  "small."  This 
will  not  be  the  situation  encountered  in  re-entry. 
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APPENDIX 


Of  the  codes  used  as  exeunples  in  this  report,  a fair 
number  are  (or  are  derived  from)  so-called  maximal  length  shift 
register  codes.  They  are  shift  register  codes  because  they  can 
be  generated  by  clocking  bits  out  of  one  end  of  a shift  register. 
The  input  to  the  shift  register  is  the  modulo-two  sum  of  a select- 
ed subset  of  the  bits  contained  in  the  previous  state  of  the  reg- 
ister. The  codes  are  maximal  length  because  this  subset  of  bits, 
the  "feedback"  bits,  has  been  chosen  so  that  the  output  sequence 
does  not  repeat  itself  until  all  possible  states  of  the  register 
other  than  all  zeros  have  occurred.  Such  a code  is  then  charac- 
terized by  a shift  register  length,  a specification  of  the  feed- 
back bits,  and  the  initial  contents  of  the  register,  or  initial 
"load."  As  an  example,  consider  the  situation  shown  in  Fig.  A-l, 
which  defines  our  conventions. 

BIT  POSITIONS 1 2 NBITS 

1110010  CODE  CHIPS 
7 CHIP  BARKER  CODE 


Pig.  A-l.  Generation  of  shift  register  code. 


I 
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The  register  is  NBITS  long,  3 in  this  case.  The  initial  load  is 
111,  these  bits  filling  the  register  from  left  to  right.  The 
register  bits  are  indexed  left  to  right  from  1 to  NBITS.  The  feed- 
back bits  are  bits  2 and  3 in  this  case.  One  mod  2 adder  is  used 
here.  It  is  equivalent  to  an  exclusive  OR  operation  on  the  two 
bits.  More  generally,  more  than  two  bits  may  be  fed  back.  The 
number  of  bits  fed  back  is  even,  and  bit  position  NBITS  is  always 
included  if  the  code  is  maximal  length.  If  the  set  of  feedback 
bits  contains  an  odd  number  of  I's, the  bit  fed  back  is  a 1.  Other- 
wise it  is  a zero.  The  first  NBITS  chips  of  the  phase  code  are 
the  initial  load  bits  in  reverse  order,  and  the  length  of  the  code 
is  2 -1  chips.  In  this  example,  the  code  produced,  or  equiva- 

lently its  complement,  is  the  7-chip  Barker  code.  If  any  of  the 

NBITS  \ 

other  possible  2 -2  starting  loads  is  used  (all  zeros  is  ex-  ( 

! 

eluded) , the  code  obtained  is  the  original  code  cycled  around.  The 
autocorrelations  of  these  other  codes  (6  in  this  example)  will  be 
different.  If  a different  set  of  feedback  bits  is  chosen,  and  the 
register  length  held  fixed,  a completely  different  set  of  codes  will 
be  generated,  not  necessarily  maximal  length. 

With  this  background.  Table  A-l  summarizes  the  codes  used 
in  this  memo  which  are  either  binary  maximal  length  shift  register 
codes  or  are  derived  from  such.  In  addition,  we  used  several  codes 
that  were  not  shift  register  derived.  The  first  of  these  was  Delong's 
31  chip  code having  lowest  peak  sidelobe  (-20.3  dB)  of  any  31  chip 
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binary  code.  The  code  bits  are 


1000011010011111101010011101110  . 


As  a check,  the  code  in  octal  form  should  be  10323752356^. 

The  other  codes  that  were  not  shift  register  derived 
were  the  two  101  chip  codes  of  Reference  2.  These  are  represented 
below,  where  each  row  is  read  left  to  right  and  rows  are  read  top 
to  bottom. 

(A)  Sequence 


+ - - + + 

+ + - + + 

+ + - + + 

- + + + - 

+ - + - + 

+ + + + + 

---  + + 
- + + + + 

+ + - - + 

•f 


- - + + + 

+ - + + - 

- + + - + 

+ 

- + - + - 

- - - + + 

+ + --- 
- - - + + 

+ + - - + 


42 


SHIFT  REGISTER  AND  RELATED  CODES 
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